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Results obtained using a recently developed calculational method for 
determining the nucleon-meson cascade induced in thick materials by high- 
energy nucleons and charged pions are presented. The calculational method 
uses the intranuclear-cascade-evaporation model to treat nonelastic collisions 
by particles with energies % 3 GeV and an extrapolation model at higher 
energies. The following configurations are considered: (a) 19.2-GeV/c 

protons incident on iron, (b) 30.3-GeV/c protons incident on iron, (c) solar 
and galactic protons incident on the moon, and (d) galactic protons incident 
on tissue. For the first three configurations, experimental results are 
available and comparisons between the experimental and calculated results 
are given. 


INTRODUCTION 

Nucleon-meson transport calculations that uti- 
lize Monte Carlo techniques in conjunction with the 
intranuclear-cascade-evaporation model for treating 
nonelastic collisions were previously carried out 
for a variety of problems involving- solar protons 
and other sources (refs. 1 to 3, for example). The 
calculational method employed for these problems is 
restricted to particle energies 3 GeV because of 
limitations imposed by the particular intranuclear- 
cascade model used (ref. ^) . To enable transport 
calculations to be performed for galactic protons 
and other high-energy sources , the calculational 
method has been extended to higher energies by us- 
ing an extrapolation model (ref. 5) to obtain the 
description of products from nonelastic collisions 
£ 3 GeV. This extrapolation method uses the dif- 
ferential cross sections for the production of nuc- 
leons and pions from 3-GeV nucleon— nucleus and 
2.5-GeV pion-nucleus collisions as predicted by the 
intranuclear-cascade-evaporation model, together 
with energy, angle, and multiplicity scaling rela- 
tions that are consistent with the sparse experi- 
mental data available for high-energy interactions, 
to estimate the particle production at the higher 
energies . 
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This recently developed high-energy transport 
method has been applied to compute the nucleon- 
meson cascade induced in thick targets for several 
source-geometry configurations, and the purpose of 
the present paper is to present some of the results 
obtained. The following cases are considered: 

(a) 19.2-GeV/c protons incident on a thick iron 
target, (b) 30.3-GeV/c protons incident on a thick 
iron target, (c) galactic and solar protons inci- 
dent on the moon, and (d) galactic protons incident 
on tissue. For the first three cases, experimental 
data exist so that comparisons with the calculated 
results can be made in order to check the validity 
of the calculational method. 

RESULTS 

19.2-GeV/c Protons Incident on Iron 

Citron et al. (ref. 6) have measured the longi- 
tudinal and lateral development of the nucleon-meson 
cascade induced in a thick iron target by a narrow 
beam of 19.2-GeV/c protons. Figure 1 shows a com- 
parison of the experimental results, the results 
from the present calculations, and the results 
available from calculations made by others for the 
longitudinal development of the cascade. The cas- 
cade development is given in terms of both the star 
density and track intensity. For the present cal- 
culations, the star density was calculated using 
the cross sections for nuclear star production in 
emulsion given by Hess et al. (ref. 10), and the 
track intensity was taken to be the flux due to pro- 
tons above 500 MeV and charged pions above 80 MeV. 
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In the experiment, the iron absorber consisted of 
slabs of iron of various thicknesses with air gaps 
between the slabs. Results for the present calcu- 
lations in figure 1 are given with and without 
these air gaps taken into account. As indicated 
in figure 1, the results of the present calcula- 
tions for the longitudinal development of the cas- 
cade agree reasonably well with the experimental 
results and with the results from other calcula- 
tions . 



FIGURE 1.— Star density and track intensity versus 
depth for 19.2-GeV/c protons incident on iron. 
Results from the present calculations are com- 
pared with the straightahead calculations of 
Alsmiller and Barish (ref. 7), the Monte Carlo 
calculations of Ranft (ref. 8), the calculations 
of O'Brien using Passow’s approximation (ref. 9), 
and the measurements of Citron et al. (ref. 6). 


The lateral distribution of the track inten- 
sity at various depths has also been calculated 
and compared with the measurements of Citron et al. 
The calculated and experimental results for the 
lateral distributions are not in good agreement. 

The reason for this discrepancy is not known at 
present. However, as discussed later, the calcu- 
lated lateral development of the cascade is in good 
agreement with other experimental data (ref. 11 ) 
for 30.3-GeV/c protons incident on iron. In prin- 
ciple, there is no reason that the calculational 
method should be less accurate for 19.2-GeV/c in- 
cident protons than for 30.3-GeV/c incident protons. 


30.3-GeV/c Protons Incident on Iron 

The lateral development of the nucleon-meson 
cascade induced in a thick iron target by a narrow 
beam of 30.3-GeV/c protons has been determined ex- 
perimentally by Awschalom et al. (ref. 11 ). The 
quantity measured was the 18 F production in alumin- 
um foils placed at various depths in the iron. 
Calculations have been carried out to determine the 
energy and spatial distributions of the nucleon and 
charged-pion fluxes for this configuration. These 
fluxes and available cross sections for producing 
18 F from aluminum (ref. l) were then used to com- 
pute the 18 F production as a function of distance 
from the beam axis. The energy dependence of these 
cross sections is such that, roughly speaking, the 
18 F production can be interpreted as a measure of 
the nucleon and charged-pion fluxes above ^ 50 MeV. 

A comparison of the calculated and experi- 
mental results is shown in figure 2 for several in- 
termediate depths in the iron. The agreement is 
quite good although the calculations give a con- 
sistently higher production near the beam axis. 
However, because of the small area and alignment of 
the aluminum foils used in the experiment, the ex- 
perimental 18 F production is expected to yield an 
underestimate of the actual production very near 

ORNL-CWG 71 -1609 


DEPTH = 
162 g/cm 2 


243 g/cm 2 


324 g/cm 2 


0 100 200 300 400 

DISTANCE FROM BEAM AXIS (g/cm 2 ) 

FIGURE 2.— Calculated and measured (ref. 11 ) 
lateral distribution of 18 F production in alu- 
minum foils at various depths in iron for 30.3- 
GeV/c incident protons. 
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the beam axis (ref. 12). The agreement at the in- 
termediate depths shown in figure 2 is fairly typi- 
cal of the agreement obtained for other comparisons 
that have been made at depths ranging from 0 to 
1045 g/cm 2 . The laterally integrated 18 F production 
is not available from the experiment, so the calcu- 
lated longitudinal development of the cascade cannot 
be compared with experiment for this case. 

Galactic and Solar Protons Incident on the Moon 

The depth-dependent radionuclide activity in- 
duced in the moon by galactic and solar proton bom- 
bardment has been calculated and compared with 
measurements made on Apollo 11 and 12 samples. As 
an example of the calculated results obtained, fig- 
ure 3 shows the depth dependence of the 26 A1 activ- 
ity. The total 26 A1 activity (i.e., the sum of the 
solar and galactic contributions) is shown in fig- 
ure 4, together with the measured activity in two 
moon rocks . Additional comparisons between calcu- 
lated and measured activities are given elsewhere 
(ref. 15). 



FIGURE 3.— 26 A1 activity versus depth in moon due 
to solar cosmic rays (SCR) and galactic cosmic 
rays (GCR) at solar minimum and solar maximum. 
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FIGURE 4.— Comparison of calculated and measured 
(refs. 13, l4) 26 A1 activity induced in the moon 
by cosmic -ray bombardment. 

The depth dependence of the neutron spectra in 
the moon has also been calculated. The thermal- 
neutron flux in the moon has been inferred from the 
Apollo measurements (ref. 16) and is in good agree- 
ment with the calculated flux. 

The same calculational method as used here to 
determine the induced activity in the moon has been 
applied recently by Gabriel and Santoro to determine 
the induced activity in the soil surrounding a 500- 
GeV proton accelerator (ref. IT). 

Galactic Protons Incident on Tissue 

The absorbed dose and dose equivalent induced 
by galactic protons in an infinite slab of tissue 
30 cm thick have been calculated using the high- 
energy transport code. An isotropic flux of galac- 
tic protons having a solar minimum energy spectrum 
(ref. 15) was taken to be incident on one side of 
the tissue slab. The nucleon-meson cascade in the 
tissue produced by incident protons in the energy 
range from 30 MeV to 200 GeV was calculated, and 
the depth dependence of the absorbed dose and dose 
equivalent was determined as described in reference 
18 . The electron-photon cascade resulting from the 
decay of neutral pions and the electrons and posi- 
trons from muon decay were taken into account in an 
approximate manner (ref. 18). The following com- 
position (in at. %) was used for the tissue: H, 

63.3; 0, 25.8; C, 9.5; N, 1.4. A tissue density of 
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1.0 g/cm 3 was used. 

The absorbed dose and average quality factor 
as a function of depth in the tissue are shown in 
figure 5. The average quality factor was computed 
by dividing the dose equivalent in a given depth 
interval by the absorbed dose in that interval. In 
determining the dose equivalent, the quality fac- 
tors for each type of particle as a function of 
linear energy transfer were obtained in the same 
manner as described in reference 17. From figure 5, 
the absorbed dose for an omnidirectional (over 2tt 
solid angle) flux incident on one side of the tis- 
sue slab is approximately h rad/y at all depths. 

For a Utt incident flux (i.e., 2tt incident flux on 
both sides of the slab), the dose would be 8 rad/y, 
which is in agreement with the value of 10 rad/y 
quoted by others (refs. 19 and 20) for a Itt flux at 
solar minimum. 
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FIGURE 5 .— Absorbed dose and average quality factor 
in tissue from galactic proton bombardment. 

The contributions to the absorbed dose and 
dose equivalent, respectively, from various kinds 
of particles are shown in figures 6 and 7« The 
total dose is obtained by adding all of the contri- 
butions shown in the figures. The histogram lab- 
eled "primary ionization” in figure 6 gives the ab- 
sorbed dose from the excitation and ionization of 
atomic electrons by those incident protons which 


have not undergone nuclear collision. The histogram 
labeled "secondary protons" gives the absorbed dose 
from the excitation and ionization of atomic elec- 
trons by protons produced from nonelastic nucleon- 
nucleus and pion-nucleus collisions and from the 
elastic collisions of nucleons and pions with hy- 
drogen nuclei. The histogram labeled "heavy nuclei" 
gives the absorbed dose from particles with mass 
number greater than one produced from nonelastic 
nucleon-nucleus and pion-nucleus collisions and the 
absorbed dose from the recoiling nuclei produced 
from elastic neutron- nucleus collisions and from 
nonelastic nucleon-nucleus and pion-nucleus colli- 
sions. The histogram labeled "charged pions" gives 
the absorbed dose from the excitation and ioniza- 
tion of atomic electrons by both positively and 
negatively charged pions produced from nucleon- 
nucleus and pion-nucleus nonelastic collisions. 

The histogram labeled "photons from neutral pions" 
gives the absorbed dose from the electron-photon 
cascade produced by the photons which arise from 
the decay of neutral pions. The histogram labeled 
"electrons, positrons, and photons" gives the 


fipwi -niA/r: 74-t«n* 



0 5 10 15 20 25 30 

DEPTH IN TISSUE (cm) 


FIGURE 6.— Contribution by particle type to the ab- 
sorbed dose in tissue from galactic proton bom- 
bardment . 
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absorbed dose from the electrons and positrons pro- 
duced by muon decay and the absorbed dose from the 
photons produced by nucleon-nucleus and pion-nucleus 
nonelastic collisions. The histogram labeled "muons" 
gives the absorbed dose from the excitation and ion- 
ization of atomic electrons by both positively and 
negatively charged muons . In figure 7 the histo- 
grams have similar meanings but give the dose equiv- 
alent from these various kinds of secondary paj*- 
ticles v 



FIGURE 7*— Contribution by particle type to the 
dose equivalent in tissue from galactic proton 
bombardment . 

It should be realized that only the protons in 
the cosmic radiation have been considered in calcu- 
lating the dose. Cosmic-ray nuclei heavier than 
protons can contribute significantly to the dose, 
especially at small depths. At present, the trans- 
port code used for these calculations is not capable 
of transporting particles heavier than nucleons. 
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